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Abstract 

A scalar potential of the form XabH^t'Pb is bounded from below if its 
matrix of quartic couplings Xab is copositive - positive on non-negative vec- 
tors. Scalar potentials of this form occur naturally for scalar dark matter 
stabilised by a Z2 symmetry. Copositivity criteria allow to derive analytic 
necessary and sufficient vacuum stability conditions for the matrix Xab- 
We review the basic properties of copositive matrices and analytic crite- 
ria for copositivity. To illustrate these, we re-derive the vacuum stability 
conditions for the inert doublet model in a simple way, and derive the 
vacuum stability conditions for the Z2 complex singlet dark matter, and 
for the model with both a complex singlet and an inert doublet invariant 
under a global U{1) symmetry. 

1 Introduction 

Interest in extended scalar sectors has been rekindled by the inconclusive results 
of a 125 GeV Higgs boson from LHC and TeVatron [1]. In the Standard Model, 
the Higgs potential is not stable up to the scale of Grand Unified Theories 
(GUT) if the Higgs boson mass is below 128 GeV, because the Higgs boson 
quartic coupling A runs to negative values below the GUT scale [2]. Extending 
the model by scalars can improve vacuum stability [3, 4, 5, G]. 

To ensure a stable vacuum, the scalar potential has to be bounded from 
below. In the Standard Model it is enough to have a positive Higgs boson quartic 
coupling A > 0. In theories with more scalar degrees of freedom, such as the 
two-Higgs-doublet model, the potential should be bounded from below in all 
directions in the field space as the field strength approaches infinity. 

Finding vacuum stability conditions is a difficult problem in general. How- 
ever, it is possible that the quartic part of the scalar potential is a quadratic 
form Xabfa'fil of squares of real fields. A Z2 symmetry that stabilises scalar dark 
matter imposes this condition in a natural way. In that case, definite analytic 
conditions for vacuum stability can be given. 

In the calculations of vacuum stability we can ignore any terms with dimen- 
sionful couplings - mass terms or soft terms -, since in the limit of large field 
values, c? < 4 terms are negligible in comparison with the quartic part of the 
scalar potential V4. Demanding that V4 > as the fields (^^ — >■ 00 is the require- 
ment of strong stability. V4 means stability in the marginal sense (then one 
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must have V2 ^0). Below we give all conditions in the marginal sense; they can 
be made strong by making the inequalities strict. 

If the quartic part of the scalar potential V4 is a biquadratic form Xab^a^l 
of real fields or gauge orbit variables, then its domain is not the whole M", 
but the non- negative orthant M" . In that case, the potential is positive if the 
matrix Xab is copositive - positive on non-negative vectors ('copositive' is short 
for 'conditionally positive'). 

The notion of copositivity was first introduced by Motzkin [7] in 1952. Hith- 
erto it has been used in the field of quadratic optimisation (see e.g. [8]). The set 
of copositive matrices is larger than the set of familiar positive (semi)definite 
matrices and includes it. While is easy to find the positive definite part of the 
parameter space using Sylvester's criterion on M", the criteria for copositivity 
are more complex. In general, testing for copositivity is NP-hard [')]. Imposing 
symmetries on the potential can considerably simplify the criteria. 

While the notion of copositivity is well known in the field of optimisation and 
nonlinear programming, to our knowledge we are the first to use these criteria 
to determine constraints on quartic couplings from vacuum stability. Recent 
reviews on copositive matrices are [10, 11] and in the Ph.D. thesis of Bundfuss 
[12]. A short review of many specific algorithms to detect copositivity is [13]. 

The purpose of this paper is to introduce basic properties of copositive ma- 
trices and criteria for copositivity. We compare copositivity with the familiar 
positive definiteness. We give explicit criteria for copositivity of 2 x 2 and 3x3 
matrices, and some general criteria in terms of the cofactors of the matrix or in 
terms of its eigenvalues and eigenvectors. We have left out lengthy proofs of the 
theorems we give and refer the interested reader to the cited papers. 

After a review we put the criteria to practice. We re-derive the vacuum 
stability conditions for the inert doublet model [14, 15] in a parameterisation 
where they are self-evident. Then we derive the necessary and sufficient vacuum 
stability conditions (25) for the full scalar potential of the Higgs boson and 
complex singlet dark matter stabilised by a Z2 symmetry. Previously, only a 
simpler potential with a global U{1) has been considered [16, 17]. We see that 
in case of singlets, in general their real components, not the polar form, have 
to be used to capture the whole parameter space. Last, we derive the vacuum 
stability conditions for a model where the scalar dark sector comprises both the 
inert doublet and a complex singlet [18, 19]. 

In Section 2, we introduce the notion of copositivity, give some criteria for 
matrices to be copositive, and compare them to the usual criteria of positive 
definiteness. In Section 3, we re-derive the vacuum stability conditions for the 
inert doublet model. In Section 4, we derive the vacuum stability conditions 
for the model of Z2 complex singlet dark matter. To give another example, we 
find vacuum stability conditions for the singlet and doublet model with a global 
U{1) symmetry in Section 5. We conclude in Section 6. 

2 Copositive Matrices 

Wc first briefly recall properties of the usual positive matrices. A symmetric 
matrix A is said to be positive semidefinite if the quadratic form Ax ^ for 
all vectors x in R" . - A positive matrix A is positive definite if the quadratic 
form Ax > for all vectors x. 
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There are several equivalent criteria for a matrix A to be positive. For ex- 
ample, all its eigenvalues have to be non-negative, as we can go to a basis where 
the matrix A is diagonal. If all Xi save one are zero, then in that basis an ^ 
for the matrix to be positive. For practical purposes, however, this is not a 
good criterion: already the eigenvalues of a 3 x 3 matrix are rather complicated 
analytical expressions. 

A more convenient one is Sylvester's criterion (see e.g. [20]). It states that 
for a symmetric matrix A to be positive semidefinite, the principal minors of A 
have to be non-negative. (The principal minors are determinants of the principal 
submatrices. The principal submatrices of A are obtained by deleting rows and 
columns of A in a symmetric way, i.e. if the ii, . . . , rows are deleted, then the 
ii, . . . ,ik columns are deleted as well. The largest principal submatrix oi A is A 
itself.) Thus if the matrix A is positive, all of its submatrices, in particular the 
diagonal elements an have to be non-negative. 

The definition of a copositive matrix is similar to the definition of a posi- 
tive definite matrix. A symmetric matrix A is copositive if the quadratic form 
x'^ Ax ^ for all vectors x ^ in the non-negative orthant M^. (Notation x ^ 
means that Xi ^ for each i = 0,...,n.) - A copositive matrix A is strictly 
copositive if the quadratic form x^ Ax > for all vectors x > 0. 

Positive matrices are a subset of copositive matrices. It is also obvious that a 
non-negative matrix A with aij ^ of any order is copositive. It is easy to show 
that for n = 2 copositive matrix is either positive-semidefinite or non-negative. 
In fact every copositive matrix of order n = 3 or n = 4 can be expressed as 
the sum + of a positive semidefinite matrix S and a nonnegative matrix TV 
[21]. For n > 4, matrices of the form S + N are a strict subset of the copositive 
matrices.^ However, because for large matrices the general tests are slow, it can 
still be useful to test if a given matrix has the form S + N [22] . 

If one component Xi is set to zero, then the quadratic form x'^Ax becomes 
a quadratic form of the remaining variables. Therefore if a matrix A of order n 
is copositive, each principal submatrix of A of order n — 1 is also copositive and 
so forth. In particular an ^ for all i. If A is strictly copositive then an > 0. 

A symmetric matrix A of order 2 is copositive if and only if [23] 



This can be easily seen: if ai2 < 0, then 011022 ^ 012 is required, but if 012 ^ 0, 
there are no constraints on an and 022 besides being positive: A is indeed 
either positive semidefinite or non-negative. Combining these statements with 
an ^ 0,a22 ^ gives (2). 

Because (2) has to hold for all 2 x 2 principal submatrices of A, we have 
a.ij + ^Jonajj ^ for all i, j. If a^i = 0, then for all j. 

Note that unlike positive definiteness, copositivity is not invariant under 
basis transformations. Consider [12] 



ai2 + \/aiia22 ^ 0. 



an 55 0,022 > 0, 



(1) 
(2) 




(3) 



^Thus random generation of copositive matrices is rather easy for n < 4, since by Cholcsky 
decomposition every positive semidefinite matrix can be factorised as S = LL^ , where L is a 
lower triangular matrix. 
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The matrix 

U^AU (4) 




is not copositivc as it docs not satisfy (2). However, copositivity is invariant 
under permutations and scaling with positive numbers. 

A symmetric matrix A of order 3 is copositivc if and only if [24, 25] 



an ^ 0,022 ^ 0,033 ^ 0, 



ai2 = ai2 + vonoii ^ 0, 



ai3 = ai3 + ^011033 ^ 0, 



(5) 



023 = a23 + V«22a33 ^ 0, 

and 

V^aiia22a33 + ai2^/a^ + 013^0^ + 023^/011 + \/2ai2ai3a23 ^ 0. (6) 
The last criterion in (6) is a simplification of 



^011022033 + ai2%/a33 + ai3\/a^+ a23\/aii ^ 0, (7) 
det A = 011022033 - (012033 + 0^3022 + 011O23) + 2012013023 > 0, (8) 

where one or the other inequality has to hold [24] . (The condition det ^ ^ is 
also part of Sylvester's criterion for positive semidefiniteness.) The conditions 
(5) simply state that the three 2x2 principal submatrices of A are copositivc. 
Only (6) looks more peculiar. 

In the literature there are also explicit criteria for matrices of order 4 and 5 
[2G] and [27] (also gives a recursive algorithm for matrices of order n), but these 
are too lengthy to reproduce here. 

For larger matrices, it is easier to use more general criteria such as the 
Cottle-Habctlcr-Lcmkc theorem [2b]: Suppose that the order n — 1 principal 
submatrices of a real symmetric matrix A of order n are copositivc. In that case 
A is non-copositive if and only if 

detA<0 and adjA^O. (9) 

The adjugate of A is the transpose of the cofactor matrix of A: 

(adj^)y = (-l)'+^Mj„ (10) 

where My is the (i, j) minor of A, the determinant of the submatrix that results 
from deleting row i and column j of A. 

To illustrate the theorem, consider the 2x2 symmetric matrix 

A.=h (11) 

\ai2 a-2,2) 

Its adjugate is 

adjA2= I I . (12) 
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By assumption, an ^ 0,022 ^ 0. The first condition det A2 = aiia22 — 0,12 
gives by 

a-b^^0^y/7i + b^0 (13) 

that if ai2 + ■\/aiiai2 < 0, then A is not copositive. Indeed, this violates (2). 
For the 2x2 matrix the second condition does not give anything new because 
ai2 may or may not be negative without violating copositivity. 

The theorem can be rephrased: Suppose that the order n — 1 principal sub- 
matrices of a real symmetric matrix A of order n are copositive. In that case A 
is copositive if and only if 

detA ^0 or some element(s) ofadj A < 0. (14) 

For A2, the condition det A2 ^ together with an ^ 0, 022 ^ is Sylvester's 
criterion of positive semidefiniteness for 2x2 matrices. Alternatively some el- 
ements of adj A2 < 0. Because an ^ 0, 022 ^ 0, it can be only —ai2 < 0. The 
condition aiia22 — 012 ^ V 012 > is exactly equivalent to (2). 

There are other algebraic criteria for copositivity. Positive definiteness can 
be checked via positivity of eigenvalues of the matrix. A similar criterion exists 
for copositivity, only now all the eigenvalues and eigenvectors of the principal 
submatrices have to be considered as well. - Kaplan's test [■_'!(] states that a 
symmetric matrix A is copositive if and only if every principal submatrix of A 
has no eigenvector v > with associated eigenvalue A < 0. (This arises from 
the positivity of the so-called Pareto eigenvalues of the matrix A [11].) 

There is also a Schur complement-like theorem for copositive matrices [30] 
stating that the matrix 

where 6 is a vector in R"^^ and C is a symmetric matrix of order n, is copositive 
if and only if 

1. a ^ 0, C is copositive, and 

2. y^{aC - bb^)y ^ for aU y £ M"-i such that b^y 0. 

The second condition is the hardest to verify. 

Because the number of principal submatrices grows exponentially with ma- 
trix size, for large matrices of order about n > 20 numerical criteria like simplex 
algorithms [31] should be used. 



3 An Example: Vacuum Stability Conditions for 
the Inert Doublet Model 

For a warmup, we re-derive the vacuum stability conditions for the inert doublet 
model [14, 15] via copositivity. The inert doublet model is a Z2-symmetric ver- 
sion of the two-Higgs-doublet model (2IIDM) [32] (see [33] for a comprehensive 
review) that contains a candidate of scalar dark matter. For completeness we 
also discuss the general 2HDM, though it is more amenable to other approaches. 
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The most general 2HDM scalar potential of two electroweak doublets Hi 
and H2 is 



(16) 



+ |i7i|2(A6ffJ-ff2 + x;hIhi) + \H2\\\7h\h2 + \;hIh{) 

= fifhl + filhj + Xihf + A2/12 + hhfhl + Xip^hlhl 
+ IA5I cos(2(/) + <j)x,)p^hlhl + 2r3r2|A6| cos((/) + (jyx,) 
+ 2rir||A7|cos(0 + 0A,), 

where we have parameterised the field bilinears as [34] 

|i?iP = /i2, 1^2!' = /»2, HIH2 = hih2pe't (17) 

The parameter |p| G [0, 1] as implied by the Cauchy inequality ^ \hIH2\ ^ 

l^l||^2|. ^ 

In the inert doublet model, the potential is required to be invariant under 
the Z2 transformations Hi — > Hi, H2 —H2, setting Ag = A7 = 0. In this case 
the quartic part of the scalar potential is a quadratic form of hi and /ij. To find 
the vacuum stability conditions, we must minimise the potential with respect to 
the free parameters p and (f). It is obvious that the minimum of the A5 term is 
given by cos(2(j) + (j)\^) = — 1. For this minimum, the matrix of quartic couplings 
for the minimal potential in the {hf, h^) basis is 

^1 i[A3 + P^(A4-|A5|)]\ 

1^1 [A3+PMA4-IA5I)] A2 )■ 

Now if A4 — IA5I ^ 0, the minimum of the potential is obtained by setting 
yO = 0, but if A4 — IA5I < 0, then the minimum is given by p = 1. Applying the 
copositivity criteria (1), (2) to (18) for these two cases yields 



Ai^O, A2^0, A3 + 2VA^^0, A3 + A4- IA5I +2yAlA^^ 0, (19) 

the well known vacuum stability conditions [14, 35, 34, 36] for the inert doublet 
model. ^ 

Derivation of the necessary and sufficient conditions (19) for the special case 
of the inert doublet model is simpler with copositivity than with the methods 
used to derive the necessary and suSicicnt positivity conditions for the most 
general 2HDM [36]. 

To apply copositivity to the latter, one would have to consider a larger 
monomial basis h= (/i^? ^1^2, /i^) which the matrix of couplings is 



A(c, 



^ Ai |A6|/9COs((?!) + 0aJ CiiA345(/3, (/>) 

I Ae I p cos ((/) + (/) As ) C2A345(p,0) |A7|pC0s((?!) + (?;>Ar) | , (20) 

\ Ci\\m5{P,4>) |A7|pCOs((?!) + (/)A7) A2 



^Had we attempted to find the vacuum stability conditions from copositivity via real com- 
ponents of the fields, we would get, besides (19), redundant inequalities like A2 + IA2I ^ 0, 
due to the SU{2) symmetry. In the basis of the eight real component fields, it would be much 
more troublesome to show that the conditions (19) are not just necessary, but also sufficient. 
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where we have denoted A345(p, c^) = A3 + [A4 + IA5I cos(2(/) + (f)\^)]. 

Now, however, the elements on the minor diagonal of the matrix (20) are 
not independent any more [37]: they depend on two parameters ci and C2 with 
the condition that 

Cl+C2 = l. (21) 

The matrices A(ci) form an affine space. 

It is easy to see that if A(cj;) is copositive for some values of Ci ~ c" satisfying 
(21) (and the whole ranges of p and (/)), then the potential (16) is bounded below. 
(Note that c° depend on the values of Ai, p and (j).) 

Namely, assuming that A(c°) is copositive, but the potential is not bounded 
below for a non- negative vector is contradictory, because the potential (16) 
is given by V{h) = h^A{ci)h for any allowed values of c,;, especially c?. Thus 
also for we have V{h-) = /i^A(c^)/i_ ^ from non-negativity of and 
definition of copositivity. 



4 Vacuum Stability Conditions for Complex Sin- 
glet Dark Matter 

A real [38, 6] or complex [16, 17, 39, 4] scalar singlet is perhaps the simplest 
candidate for dark matter. The scalar sector of the latter model comprises the 
Standard Model Higgs doublet Hi and the complex singlet 

g^SH+^^ (22) 

The most general scalar potential invariant under the Z2 transformations 
Hi Hi, S ^ -S is 

V = ^,l\Hi\' + mII^I^ + ^[S' + iS^r] + Xi\Hi\' 
+ ^s\Sf + ^[S' + {S^r] + ^\S\^ [S' + (5t)2] (23) 
+ Xsi\S\'\Hi\' + ^\Hi\^S' + {S^r]. 

Only the Hermitian square jffi |^ of the Higgs field appears since it is the only 
scalar doublet in the model, and only the squares S]j and 5^ of the components 
of the singlet appear because of the Z2 symmetry. 

The matrix of quartic couplings A in the {hf, Sj^, S*^) basis is given by 



4A 



^ 4Ai Asi + A^^ Asi — A^^ 

Asi + A'si As + X's + A'i As - SX's | • (24) 
\Asi — Agj As — 3A^ As + A'g — Ag^ 



Applying the criteria for copositivity (5) and (6) to the matrix of quartic 
couplings (24) yields 
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Ai ^5 0, Xs + X's + A'^ ^0, As + A's - A'^ ^ 0, 



Am = Asi + X'si + 2^ Xi{Xs + X's + \'^) ^ 0, 
~XiA = Asi - A'si + 2y^Ai(As + A^-A'^) ^ 0, 
~Xha = As - 3A's + y^(As + A^)2 - A'^^ ^ q, 



(25) 



+ Asi (^^As + A'g - A'^ + ^As + A'g + A'< 

+ Asi As + Ag - Ag - y^As + Ag + Ac 

+ \/2\/AiijAiAAffA ^ 0. 

Wc emphasise that these are the necessary and sufficient conditions for vacuum 
stabiUty.'' 

Previously [17], the vacuum stabihty conditions have been derived only in 
the case of an extra global J7(l) symmetry imposed on the potential (23). Then 
As = As = A's]^ = 0, and the conditions (25) reduce to 



Ai ^0, As ^ 0, Asi ^ -2v/aTa^ (26) 

that - taking into account normalisation conventions - coincide with the condi- 
tions given in [17]. In that case the last inequality in (25) becomes redundant, 
being identically true by virtue of (26). 

Could one derive the conditions (25) by parameterising the singlet field as 
S — se"^^ and considering the copositivity criteria (1) and (2) for the 2x2 
matrix 

^^i Ai i(Asi+A'siCos20s) \ ^27) 

\^i(Asi + A;^iCos20s) As + A'5Cos40s + A'^cos20sy 

of couplings? 

In the C/(l) case we just encountered, the terms dependent on 0s are for- 
bidden and (27) gives exactly the conditions (26). 

Otherwise, we have to minimise the potential with respect to 0s- Already 
for the self-couplings of S*, however, this approach fails to capture the whole 
allowed region. If h\ = 0, we must have 

As + A's cos 40s -t- A's cos 20s 0. (28) 

If A's = 0, then the minimum of As + A^ cos 40s is obviously As — |As|. Else, 
the extremum condition is 

2A'5 sin 40s -I- A's sin 20s = (A's + 4A's cos 20s) sin 20s = 0, (29) 



■^The previously given conditions in [ID] pertaining to H\ and S are necessary conditions 
for positivity, not copositivity. 
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Figure 1: Allowed parameter spaee for the quartic self-couplings of S. Area 
with dashed border: the parameter space allowed for A'^ = 0. Only X'g ^ 
is allowed by positivity (light red), while copositivity also includes the area 
below zero (light green). Area with solid border: the parameter space allowed 
by positivity (red), the constraint from minimising the potential (red and blue) 
and copositivity for real fields (the whole area) for A'^ = 1/2. 



± arccos ( —437- I + 2n7r 



giving 05 = ±n§ and 0s = 5 
reproduces A5 + A^ ± A'^ ^0. The latter solution gives 



The former solution 



\/l2 

As - A;5 - ^ ^ 0. (30) 

From Sylvester's criterion for the 23 submatrix of (24) one gets a similar 
condition for the usual positivity of singlet self-couplings: 

8(As - X's)X's - Af ^ 0. (31) 

The regions are illustrated in Figure 1. If A'^ = (area with dashed border), 
positivity only allows A5 ^ X'g ^ (light red), while copositivity allows As — 
IA5I > (light red and light green). 

For Ag ~ 1/2, the positivity condition (31) allows only the red region, the 
copositivity condition (30) with minimising the potential allows the red and 
blue regions, but a sliver (green) of the full allowed region given by A^^ ^ 
(25) is left out. The full allowed region (solid green border) encompasses the red 
and blue regions. 

We see that because there are several phase-dependent terms in the potential, 
the complex singlet has to be decomposed into real and imaginary parts instead 
of the polar form to guarantee that the obtained vacuum stability conditions are 
necessary and sufficient. An alternative would be to scan over 0s numerically. 



5 Complex Singlet and Inert Doublet with a 
Global U{1) Symmetry 

As another example we consider extending the Standard Model with both the 
complex singlet S and the inert doublet H2 [18, 19]. This is the minimal exten- 
sion of the Standard Model that enables soft cubic couplings that can induce 
electroweak symmetry breaking [40], or semi-annihilation [41]. 



The necessary and sufficient vacuum stability conditions of the general Z2- 
symmetric model [19] can be derived via the Cottle-Habetler-Lemke theorem 
from the matrix of quartic couplings in the basis {hi, Sj^, 5^).'* However, to 
present a useful illustration of the copositivity criteria, we impose a global U{1) 
symmetry on the potential. With the given field content, the scalar potential 
invariant under the global U{1) symmetry Hi Hi, H2 — >■ e^°'H2, S — )• e^^S 
transformations is 

V = ^,j\Hi\' + ^il\H2\'- + fil\S\' + ^{HlHiS + hIh^S^) 

+ Asl^l^ + \i\Hi\^ + A2|i72|' + \3\Hi\^\H2\^ + Xi{HlH2){HlHi) (32) 
+ Asi|5|2|i/i|2 + As2|5|2|i72|'. 

In the paramctcrisation (17), the quartic part of the potential (32) is 

Va = Xgs^ + \ih\ + X^hl + X^hlhl + Xip^hlhl + Xsihls^ + Xs2hls^, (33) 

where we have used S = se"^^. The potential is minimised for p = if A4 ^ 0, 
and for p = 1 if A4 < 0. As a shorthand for this, we can use the Heavyside 
theta function 0(— A4). The matrix of couplings A for the minimal potential in 
the {hf, h\, s^) basis is given by 



2A 



/ 2Ai A3 + 0(-A4)A4 A51 

A3 + 0(-A4)A4 2A2 A52 I . (34) 

\ Xsi Xs2 2A5y 



Copositivity criteria (5) and (6) yield the necessary and sufficient vacuum 
stability conditions: 



Ai ^0, A2 ^ 0, A5 ^ 0, 
A12 = A3 + 6i(-A4)A4 + 2VaIa^ ^ 0, 



Ais = Asi + 2\/XiXs ^ 0, 



(35) 



A2S = A52 + 2-\/ A2A5 ^ 0, 

and 



\/AiA2As + [A3 + 0(— A4) A4] \/As + Asi \/A2 + As2 vAi 

/ - - - (36) 
+ V A12A1SA2S ^ 0. 

Note that it is easy to take into account some terms that break the U{\) 
symmetry. For example, we can add the A'g term to the potential. As discussed 
in Section 4, this corresponds to A5 — > As ~ [A^l in (35) and (36). In the same 
way, we can add the A5 term, with A3 + (?(— A4)A4 — >■ A3 + 6'(— A4 + [AsD (A4 — IA5I) 
in (35) and (36). 



*The conditions given in [ID] arc necessary conditions for positivity, not copositivity; in [3] 
sufficient conditions were used. 
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6 Conclusions 



For a scalar potential whose d ~ A part is a biquadratie form of the fields, 
Xab'Pa'Pb^ the vacuum stability conditions for the potential to be bounded be- 
low coincide with the criteria for the matrix of its quartic couplings Xab to be 
copositive or positive on non-negative vectors. The notion of copositivity has 
found use in convex optimization, but has hitherto been overlooked in particle 
physics. Copositivity allows for a larger parameter space than positivity. For the 
latter, Sylvester's criterion gives a simple set of inequalities, but conditions for 
the former are more complicated. 

We reviewed basic properties of copositive matrices and explicit criteria for 
copositivity of 2 x 2 (these are already well known in particle physics literature, 
though not under that name) and of 3 x 3 matrices. The Cottle-Habetler-Lemke 
theorem or other algebraic criteria can be used for larger matrices. In code, the 
theorem is easier to use even for small matrices, as it is a polynomial criterion 
that does not involve calculation of square roots. For really large matrices (of 
about order n > 20), approximate algorithms have to be used. 

We illustrated this method by deriving the vacuum stability conditions for 
some dark matter models where the dark sector comprises either a complex sin- 
glet, an inert doublet, or both. In conclusion, copositivity allows to find analytic 
necessary and sufficient vacuum stability conditions for complicated models in 
an important special case. 
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